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ABSTRACT: Inspired by quantum cosmology, in which the wave function of the uni- 
verse is annihilated by the total Hamiltonian, we consider the internal dynamics of a simple 
particle system in an energy eigenstate. Such a system does not possess a uniquely defined 
time parameter and all physical questions about it must be posed without reference to 
time. We consider in particular the question, what is the probability that the system's 
trajectory passes through a set of regions of configuration space without reference to time? 
We first consider the classical case, where the answer has a variety of forms in terms of a 
phase space probability distribution function. We then consider the quantum case, and we 
analyze this question using the decoherent histories approach to quantum theory, adapted 
to questions which do not involve time. When the histories are decoherent, the prob- 
abilities approximately coincide with the classical case, with the phase space probability 
distribution replaced by the Wigner function of the quantum state. For some initial states, 
decoherence requires an environment, and we compute the required influence functional 
and examine some of its properties. Special attention is given to the inner product used 
in the construction (the induced or Rieffel inner product), the construction of class op- 
erators describing the histories, and the extent to which reparametrization invariance is 
respected. Our results indicate that simple systems without an explicit time parameter 



may be quantized using the decoherent histories approach, and the expected classical limit 
extracted. The results support, for simple models, the usual heuristic proposals for the 
probability distribution function associated with a semiclassical wave function satisfying 
the Wheeler-DeWitt equation. 
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1. INTRODUCTION 



There are a variety of interesting physical situations in which our classical view of 
the world inspires us to ask questions that quantum theory does not easily answer. For 
example, classical mechanics concerns simultaneously specified values of coordinates and 
momenta while quantum theory has to go through some contortions to say what this means 
operationally. Another important class of problems of this type are those that involve time 
in a non-trivial way. The arrival time problem and tunneling time problem, for example, 
have been the subject of considerable recent interest [1] . 

Perhaps more intriguing than these are questions in quantum mechanics that do not 
involve time at all. Consider for example the following situation. Suppose we have a 
system of particles in a state of fixed total energy. It could, for example, be a light particle 
orbiting a massive particle. Then classically, we can ask whether the light particle passes 
through a certain region of configuration space at any stage during its orbit. Or we can 
ask which of the possible classical orbits the light particle follows. The important point 
here is that these sorts of questions do not involve time explicitly. Much experimental 
and observation data is in fact of this type. For example, astronomical observations yield 
planetary orbits, and particle physics experiments often yield a photograph of a track in a 
bubble chamber. We would like to ask the same question in quantum theory: given that 
the system is in an energy eigenstate, what is the probability that it is found in a certain 
region of configuration space, irrespective of time? 

The primary motive for considering this question is quantum cosmology [2] . There, in 
simple cosmological models, the wave function of the system obeys the Wheeler-DeWitt 
equation. 



where H is the total Hamiltonian of the gravitational field plus all matter sources in the 




(1.1) 
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universe. The most significant feature of this equation is that it contains absolutely no 
reference to time whatsoever. It is usually argued that "time", or more precisely, the 
physical systems that we use to measure time, are contained already in the gravitational 
and matter fields [3,4,5]. While this is very plausible it leaves us with the question as to 
how to extract interesting physical predictions from this wave function, given the absence 
of the time coordinate that plays such a central role in standard quantum theory. In 
quantum cosmology, heuristic methods (mainly the "WKB interpretation") have been 
used to extract useful predictions [6], but what we are concerned with here is how such 
heuristic ideas may be incorporated in a properly defined interpretational framework for 
the quantum theory of timeless models. Furthermore, although simple (minisuperspace) 
models are unlikely to be reasonable physical approximations to a full quantum gravity 
theory, it does seem likely that a quantum theory of gravity will centre around a timeless 
equation of the form (1.1) (the loop variable approach, for example, involves such an 
equation [7]), hence it remains important to understand the quantization and interpretation 
of such systems at an elementary level. 

Quantum cosmology, then, is our main motivation for studying the internal dynamics 
of a closed system in an energy eigenstate. The question we shall focus on is this: given an 
energy eigenstate, what is the probability that the system enters a region A of configuration 
space? Or similarly, suppose the coordinates of the system are x = {xi,X2,-- -Xn), what 
is the probability distribution of, say, ■ ■ -Xn, given the value of xi"? Classically, these 
questions may be answered reasonably easily. We look for all the phase space initial data 
points whose classical trajectories pass through A, and then the desired probability for 
entering A is the probability measure on this subset of phase space. In quantum theory the 
question is considerably more complicated, and, like the question of phase space samplings, 
a variety of different (generally inequivalent) approaches may be employed. 

It is perhaps of interest to spell out in more detail how the absence of a time parameter 
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affects these considerations. In a standard non-relativistic particle theory, the variables of 
interest are, for example, positions at a fixed moment of time x(t). Here, the time t is 
regarded as an observable physical parameter. In a theory without time, by contrast, the 
quantities of interest are curves in configuration space x(s) (or more generally, in phase 
space). Here, s is not a physically measurable time, but is simply a parameter labeling 
the points along the curve, and the curves are parameterized in this way for mathematical 
convenience. Furthermore, one of the characteristic features of genuinely timeless theories 
is that none of the components of x are monotonic in s. This means that it is not possible 
to use one of the components of x as a "time" parameter, except for local sections of 
the curve. Despite these features, such classical theories are well-defined and predictive. 
Indeed, many classical cosmological models are of this type. For example, in the massive 
scalar field cosmological model (for a positive curvature FRW metric) [8], the classical 
solutions go backwards and forwards in both the scale factor a and the scalar field (p. 
Therefore, the most general such classical theory is one in which there is a probability 
distribution on the set of classical trajectories, and this is also well-defined (subject to 
careful normalization) as we shall see. Such probability distributions are of particular 
interest for predicting, for example, the likelihood of the initial conditions for inflation [9] . 

A closely related issue is the fact that a vanishing Hamiltonian (or in the quantum 
theory, Eq.(l.l)) is associated with the symmetry of reparametrization invariance, which 
is essentially the freedom to redefine the parameter s labeling points along a trajectory. 
Individual phase space points are not reparametrization- invariant, since they are moved 
along the classical trajectories by a reparametrization. But a useful invariant quantity is 
the entire classical trajectory, as we shall see later in more detail. This simple observation 
turns out to be a useful focal point for what we do in this paper. 

Turning now to the quantum theory, there are, as stated, many possible approaches 
to the these systems, all of which are complicated by the absence of a time parameter. 
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However, the aspects of the classical theory outlined above suggest that a particularly- 
useful approach to quantization is the decoherent histories approach [10,11,12]. This is 
because it deals directly with entire trajectories and does not obviously require a time 
coordinate. The aim of this paper is to show that the decoherent histories approach can 
be used to calculate probabilities for histories in configuration space in simple timeless 
models. In particular, we shall show that in the classical limit, the decoherent histories 
approach produces a set of classical trajectories with a probability measure on that set. 

The other significant approach to this problem that is currently being pursued involves 
constructing operators corresponding to the questions posed above, and which commute 
with the Hamiltonian, so are "observables" [13,14,15,16,17,18]. For a free particle in two 
dimensions, for example, the classical trajectories have the form, 

xiit) = xi + ^, x2{t)=X2 + — (1.2) 
m m 

and we may eliminate t between them to write, 

xi{t)^xi + ^{x2{t)-X2) (1.3) 
P2 

This is the classical answer to the question, what is the value of xi at a given value of X2? 
It commutes with the free particle Hamiltonian, 

H=l{pl+Pi) (1.4) 

so is an observable and we may find states that are eigenstates of both (1.3) and (1.4), from 
which one may begin to address the questions set out above. We will not say much about 
this approach here, but it is important to mention because any other approach ought to 
make some kind of contact with it at some stage, and in particular, the decoherent histories 
approach must contain some notion corresponding to the notion of an observable in the 
operator approach. (See also Refs. [19,20] for some much earlier approaches to these issues). 
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We begin Section 2 by describing the classical result. We introduce a classical phase 
space distribution function w(p,x) and compute the probability that a trajectory in con- 
figuration space passes through a region A. Most of the key ideas for this paper are in fact 
contained in this classical result. In particular, we discuss the reparametrization invariance 
of the system, and introduce observables corresponding to entire classical trajectories. We 
write the classical result in a number of different forms, including a form in terms of a 
flux across a hypersurface, closely related to the heuristic WKB interpretation of quantum 
cosmology. 

We begin the quantum case in Section 3, with the construction of the decoherence 
functional for timeless models. This has two important aspects. The first concerns the 
choice of inner product in the construction, since solutions to equations of the form Eq.(l.l) 
are typically not normalizable in the simple Schrodinger inner product. The appropriate 
choice is the induced (or Rieffel) inner product, which we describe. We also show how this 
inner product for the quantum case implies a useful normalization for the classical phase 
space distribution function. The second aspect is the construction of the class operators, 
which, in this case, are propagators describing coarse-grained sets of histories passing 
through restricted regions of configurations space. 

In Section 4, we discuss the semiclassical limit of the decoherence functional. A key 
step is the construction of class operators corresponding to restricted sets of histories en- 
tering the region A. The obvious candidates for these class operators are not in fact 
compatible with the constraint equation, and we therefore show how they may be appro- 
priately modified. This turns out in fact to be the crucial step in the construction of the 
decoherence functional. We then show that, for the special initial states for which the 
histories are decoherent, the probabilities for the histories approximately coincide with the 
classical case, with the phase space distribution function u'(p,x) replaced by the Wigner 
function W{p, x) of the quantum system. 
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In Section 5, we consider the special case in which the system is a collection of har- 
monic oscillators in a fixed energy eigenstate. For this system it is possible to introduce a 
special class of eigenstates of the Hamiltonian sometimes called "timeless coherent states" , 
which have the property that they are concentrated about an entire classical phase space 
trajectory. We discuss the decoherence and probabilities associated with these states and 
obtain the intuitively expected physical results for the probabilities of entering a region A. 

Since decoherence is only obtained for special initial states, we consider, in Section 6, 
the addition of an environment to produce decoherence for a wide variety of initial states. 
We repeat the calculation of decoherence and probabilities with intuitively expected results, 
in agreement with classical expectations. 

The calculations of Sections 4 and 6 used initial states consisting of single WKB wave 
packets. In Section 7, we therefore extend to the case of superpositions of such wave 
packets. This turns out in fact to be straightforward, and very similar to earlier calculations 
performed with the reduced density matrix. We easily find that the interference terms 
between different WKB wave packets are very small. We summarize and conclude in 
Section 8. 

Finally, we outline the relationship of the present paper to other works in the field. 
This paper is part of a programme to apply the decoherent histories approach to quantum 
cosmology, or models of the type used in quantum cosmology, and ultimately to quantum 
gravity more generally. The decoherent histories approach generally and its applications to 
quantum cosmology have been set out at length by Hartle in his 1992 Les Houches Lectures 
[21], and more recent relevant aspects of this were discussed by Hartle and Marolf [22]. 
This formalism have been applied to particular models in three places, most recently by the 
present authors [23] , who used it to construct in detail the decoherence functional for the 
relativistic particle, and compute probabilities for crossing spacelike surfaces. The present 
paper is in some ways an extension of that work. Whelan [24] has used the formalism 
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to compute probabilities on timelike surfaces for the relativistic particle. Also, Craig and 
Hartle [25] have applied the formalism to a Bianchi IX quantum cosmological model. The 
last two papers use the Klein-Gordon inner product whereas here, we use the positive- 
definite induced inner product to construct the decoherence functional. There is also some 
connection with the work on probabilities for non-trivial spacetime coarse grainings in 
non-relativistic quantum mechanics [26]. 

It is also perhaps worth mentioning that this paper is very much in the spirit of Ref.[27], 
which attempts to interpret the Wheeler-DeWitt equation in terms of emergent trajecto- 
ries, by introducing model detectors into the Hamiltonian. This was in turn inspired by 
Barbour's observations [28], on the similarity between the Wheeler-DeWitt equation and 
Mott's calculation showing the emergence of a straight line track from a spherical wave 
in alpha decay [29] , together with some of Barbour's more general observations about the 
Wheeler-DeWitt equation and timeless theories [28,30,31]. 



2. THE CLASSICAL CASE 

We are interested in the question, "What is the probability associated with a given 
region of configuration space when the system is in an energy eigenstate?" . We begin by 
analyzing the classical problem. 

We will consider a classical system described by a 2n-dimensional phase space, with 
coordinates and momenta (x, p) = (xj^jPf^), and Hamiltonian 

« = E l^ + i'w (2-1) 

k=l \ / 

More generally, we are interested in a system for which the kinetic part of the Hamiltonian 
has the form g^^{^)PkPji where g^^{^) is an inverse metric of hyperbolic signature. Most 
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minisuperspace models in quantum cosmology have a Hamiltonian of this form. However, 
the focus of this paper is the timelessness of the system, and the form of the configuration 
space metric turns out to be unimportant. So for simplicity, we will concentrate on the 
form Eq.(2.1). 

We assume that there is a classical phase space distribution function w{p,x.), which is 
normalized according to 

J dy dJ^x w(p, x) = 1 (2.2) 

and obeys the evolution equation 

where { , } denotes the Poisson bracket. The interesting case is that in which w is the 
classical analogue of an energy eigenstate, in which case dw/dt = 0, so the evolution 
equation is simply 

{H,w} = (2.4) 

It follows that 

«;(p^'(t),x^^(t)) = «;(p(0),x(0)) (2.5) 

where p'^'(t), x^^(t) are the classical solutions with initial data p(0),x(0), so w is constant 
along the classical orbits. (The normalization of w then becomes an issue if the classical 
orbits are infinite, but we will return to this in the quantum case discussed below). 

Given a set of classical solutions (p'^'(t), x'^'(t)), and a phase space distribution function 
w, we are interested in the probability that a classical solution will pass through a region 
A of configuration space. We construct this as follows. First of all we introduce the 
characteristic function of the region A, 

/a(x) = |1' '{^^^'^ (2.6) 
•'^^ ^ lO otherwise. ^ ^ 
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To see whether the classical trajectory x^'(t) intersects this region, consider the phase 
space function 

/oo 
dt5H(x-x^^(t)) (2.7) 
-00 

(In the case of periodic classical orbits, the range of t is taken to be equal to the period). 
This function is positive for points x on the classical trajectory labeled by pg, xq and zero 
otherwise. Hence intersection of the classical trajectory with the region A means, 

[ cZ^x /a(x) r dt 5(")(x - x^^(t)) > (2.8) 

Or equivalently, that 

/GO 
dtf^{^'\t))>0 (2.9) 
-00 

This quantity is essentially the amount of parameter time the trajectory spends in the 
region A. We may now write down the probability for a classical trajectory entering the 
region A. It is, 

PA = j (f'pod^xo «;(po, xo) e (^J^^ dt fAi^'\t)) - (2.10) 

In this construction, e is a small positive number that is eventually sent to zero, and is 
included to avoid possible ambiguities in the 6'-function at zero argument. The 6*- function 
ensures that the phase space integral is over all initial data whose corresponding classical 
trajectories spend a time greater than e in the region A. 

The classical solution x'^'(t) depends on some fiducial initial coordinates and momenta, 
Xq and po, say. In the case of a free particle, for example, 

x^^(t) = xo + ^ (2.11) 

The construction is independent of the choice of fiducial initial points. If we shift xg, 
Po along the classical trajectories, the measure, phase space distribution function w and 
the ^-function are all invariant. Hence the integral over xq, po is effectively a sum over 
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classical trajectories. The shift along the classical trajectories may also be thought of 
as a reparametrization, and the quantity (2.10) is in fact a reparametrization-invariant 
expression of the notion of a classical trajectory. This means that the probability (2.10) 
has the form of a phase space overlap of the "state" with a reparametrization-invariant 
operator. 

It is useful also to write this result in a different form, which will be more relevant to 
the results we get in the quantum theory case. In the quantum theory, we generally deal 
with propagation between fixed points in configuration space, rather than with phase space 
point. Therefore, in the free particle case, consider the change of variables from xo,po to 
xq, Xj^, where 

^/ = ^0 + j^r (2.12) 

Hence xy is the position after evolution for starting from xq for parameter time r. The 
probability then becomes 

PA = ^J d^xfd'^xo «;(po,xo) d (^J^ dt /A(x^(t)) - (2.13) 

where po = M(xy — xo)/t and 

4it) = ^0+^^^^^^t (2.14) 

The parameter r may in fact be scaled out of the whole expression, hence the probability 
is independent of it. 

The result now has the form of an integral over "initial" and "final" points, analogous 
to similar results in quantum theory. The result is again essentially a sum over classical 
trajectories with the trajectories now labeled by any pair of points xq, xy along the trajec- 
tories, and is invariant under shifting xq or xy along those trajectories. Naively, one might 
have thought that the restriction to paths that pass through A is imposed by summing 
over all finite length classical paths which intersect A as they go from the "initial" point 
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xq to "final" point xj-, that is, A lies between the initial and final points. This is also 



what one might naively expect in the quantum theory version. However, one can see from 
the above construction that the correct answer is in fact to sum over all classical paths 



along the entire trajectory, even if A does not lie between the two points (see Fig.l). This 
feature is related to the reparametrization invariance of the system. 

The above point turns out to be quite crucial to what follows in the rest of this paper, 
so it is worth saying it in an alternative form. Loosely speaking, the statement is that 
only the entire classical path respects the reparametrization invariance associated with the 
constraint equation. A section of the classical path does not. This may be expressed more 
precisely in terms of the function ^(x, po,xo) introduced in Eq.(2.7). This function is 
concentrated on the entire classical trajectory, and is zero when x is not on the trajectory. 
It is easy to see that it has vanishing Poisson bracket with the Hamiltonian H = i?(po, xg), 
since we have 



This is the precise sense in which the entire trajectory is reparametrization invariant, 
and the phase space function A may be regarded as an observable - a quantity which 
commutes with the constraint H [16,17]. By way of comparison, consider a second phase 
space function similarly defined, but on only a finite section of trajectory. 



(which can be of infinite length) passing through xq and xy that intersect A at any point 




= 



(2.15) 




(2.16) 



It is easily seen that 



{ii',S(x, po,xo)} 



(5(x-x^X^)) + <^(x-x^'(0)) 



(2.17) 



13 



Hence B "almost" commutes with H, failing only at the end points, and it is in this sense 
that a finite section of trajectory does not fully respect reparametrization invariance. 

A third version of the classical result is also useful. It is of interest to obtain an 
expression for the probability for intersecting an {n — 1) -dimensional surface S. Since the 
result (2.10) involves the parameter time spent in a finite volume region A it does not apply 
immediately. However, suppose that the set of trajectories contained in the probability 
distribution w intersect the (n — l)-dimensional surface E only once. Then we may consider 
a finite volume region A obtained by thickening S along the direction of the classical fiow. 
If this thickening is by a small (positive) parameter time At, then the quantity appearing 
in the ^-function in (2.10) is 

J dtj^(^^ 5(")(x-x^^(t)) = At J dt^d^-^x n-^^ 5(")(x-x'=^(t)) 

= At/[S,x^'(t)] (2.18) 

where n is the normal to E, and we suppose that the normal is chosen so that n • dx.^^/dt 
is positive. The quantity /[S, x'^^(t)], in a more general context, is the intersection number 
of the curve x'^^(t) with the surface S, and takes the value for no intersections, or ±1 
(depending on whether there is an even or odd number of intersections). In this case we 
have assumed that the trajectories intersect at most once, hence / = or 1. We then have 

e{AtI - e) ^ e{I - e') ^ I (2.19) 

(where e = Ate') and the probability for intersecting S may be written 

Pj:=JdtJ rpod'^xo «;(po,xo) j^d^'^xn-^^^ 5(")(x - x^^(t)) (2.20) 

At each t, we may perform a change of variables from p,x to new variables p' = P^\t), 
x' = x'^^(t), and using Eq.(2.5), we obtain the result 

P^ = Jij^^j^ c^V d^-\' n • p' «;(p', x') (2.21) 
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Finally, the integrand is now in fact independent of t, so the t integral leads to an overall 
factor. (This might be infinite but is regularized as discussed below). We therefore drop 
the t integral. 

This result is relevant for the following reason. In the heuristic "WKB interpretation" 
of quantum cosmology, one considers WKB solutions to the Wheeler-DeWitt equation of 
the form 



It is usually asserted that this corresponds to a set of classical trajectories with momentum 
p = V^", and with a probability of intersecting a surface E given in terms of the fliix of 
the wave function across the surface [6,9]. As we shall show, from the decoherent histories 
analysis, the quantum theory gives a probability for crossing a surface E proportional to 
Eq.(2.21) with w replaced by the Wigner function of the quantum theory. The Wigner 
function of the WKB wave function is, approximately [32], 



Inserting in Eq.(2.21), we therefore obtain, up to overall factors, the probability distribu- 
tion. 



* = Ce 



.iS 



(2.22) 



W-(p,x) = |C(x)|2 5(p-V5) 



(2.23) 




(2.24) 



We therefore have agreement with the usual heuristic analysis. 



3. THE QUANTUM CASE 



3(A). Decoherent Histories Approach 



The decoherent histories approach to quantum theory is described at length elsewhere 
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[10,11,12], so only the briefest review will be given here. The central object of interest is 
the decoherence functional, 

D{a,a!) = Tr(CapCl,') (3.1) 
where the histories are characterized by the class operators C^, which satisfy 

^Ca = l (3.2) 

a 

and therefore 

J^Dia.a) = Trp=l (3.3) 

a,a' 

In non-relativistic quantum mechanics, the class operators are given by time-ordered se- 
quences of projection operators 

Ca = PaM---Pa^{tl) (3.4) 

(and by sums of terms of this form), where a denotes the string of alternatives ai,a2 - ■ - an- 
The theory is however more general than this and we will exploit this generality here. 

Intuitively, the decoherence functional is a measure of the interference between pairs 
of histories a, a'. When its real part is zero for a ^ a', we say that the histories are 
consistent and probabilities 

p{a) = D{a,a) (3.5) 

obeying the usual probability sum rules may be assigned to them. Typical physical mech- 
anisms which produce this situation usually cause both the real and imaginary part of 
Z)(a, a') to vanish. This condition is usually called decoherence of histories, and is related 
to the existence of so-called generalized records [10,33]. 

In the non-relativistic case, for histories characterized by projections onto configuration 
space, a path integral version of the decoherence functional is available, and can be very 
useful. It has the form, 

D{a,a') = jj^^ //^^""""P (^'^[''^^^^ ~ ^"^[^^^^0 '"^'"^'^^^ ^^-^^ 
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where the sum is over pairs of paths x{t), y{t) passing through the pairs of regions a, a!. 
This is equivalent to the form (3.1), (3.4), when the histories are strings of projections onto 
ranges of positions. Eq.(3.6) is a useful starting point for the generalization to timeless 
theories. 

We would like to construct the decoherence functional for the situation in question, 
in which we have a system in an energy eigenstate, and we ask questions which do not 
refer in any way to time. Two new issues arise in this case. The first concerns the inner 
product through which the various piece of the decoherence functional are put together. 
The second concerns the construction of the class operators. We take each in turn. 



3(B). The Induced Inner Product 



For many situations, and especially for the analogous situation in quantum cosmology, 
the Hamiltonian has a continuous spectrum so the energy eigenstates are not normalizable 
in the usual inner product, 

(*l|*2) = j d'^x *|(x)*2(x) (3.7) 

A way to deal with this has been developed, and goes by the name of the induced inner 
product, or Rieffel induction [34,22]. Consider the eigenvalue equation 

H\^Ex)^E\^Ex) (3.8) 

where A denotes the degeneracy. These eigenstates will typically satisfy 

{^E'X'\'^EX) = S{E-E')5xy (3.9) 

from which it is clear that the inner product diverges when E = E' . The induced inner 
product on a set of eigenstates of fixed E is defined, loosely speaking, by discarding the 
(5-function 5{E — E'). That is, the induced or physical inner product is then defined by 

{'^EX'\'^Ex)phys = hx' (3.10) 
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This procedure can be defined quite rigorously, and has been discussed at some length 
in Refs. [34,22]. We will use it here to construct the decoherence functional. A simple 
prescription for using it in the decoherence functional is to regularize each propagator and 
energy eigenstate by using a different energy for each. The final answer will then involve 
a number of 5-functions in energy, as in (3.9), which are simply dropped. 

The induced inner product normalization for the wave functions does in fact suggest 
a normalization scheme for the corresponding classical phase space distribution function, 
which, recall, is not normalizable in the case where the classical trajectories are infinite 
(since w is constant along those trajectories). The idea is to consider the normalization of 
the Wigner function in the quantum case. This is defined by 



(See Refs. [35,36] for properties of the Wigner function.) For an energy eigenstate l^^^;) we 
first of all construct a regularized density operator. 




(3.11) 



with inverse 




(3.12) 



PEE' = 



(3.13) 



which is normalized by 



Tv{pEE') = S{E-E') 



(3.14) 



hence the corresponding Wigner function is normalized according to 




(3.15) 



Now notice that the density operator obeys the equation 



[H,pEE'] = {E-E')pEE' 



(3.16) 
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Taking the Wigner transform of this equation we obtain 

CWee' = i{E - E')Wee' (3.17) 

where C is the phase space operator 

f. _-sr^ ( Pk 9 dV d \ /oic\ 
^ = l.[-Mg^^+g^^g^J+^^ (3-18) 

It is a sum of the classical Liouville operator, plus a term Cq describing quantum modifi- 
cations. 

We may now see how to normalize the classical case. We take the classical distribution 
function to be described by Eq.(3.17) with the quantum term Cq set to zero. The Liouville 
operator may then be written, 

C^-l (3.19) 
for some parameter s, and Eq.(3.17) may be solved, with the result, 

WEE' = e-''^^-^'^WEE (3.20) 

The exponential factor now effectively regularizes the phase space distribution function. 
Wee is constant along the classical trajectories, but wee' is not, and, in the normalization 
(3.15), the part of the integral along the trajectories is an integral over s which produces 
the 5-function. 



3(C). Construction of the Decoherence Functional. 



The next part of the construction of the decoherence functional is the class operators, 
Ca- These are to describe histories of fixed energy which do or do not pass through the 
region A without regard to time, denoted a — A and a = A respectively. We follow 
Refs. [21,22], in part. 
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Consider first the ampfitude to go from xq at time t = to at time t = r passing 
through the region A, or not, at any time in between. This is given by 

gai^f, r|xo, 0) = ^ Vx{t) exp (^^S[x{t)]^ (3.21) 

where the sum is over all paths x{t) with t in the range [0, r] which pass through A, or 
never pass through A. It therefore satisfies, 

^9a = 9A+ 9a = 9 (3-22) 



a 



where g = g{^f, t|xo, 0) is the unrestricted propagator. There are many ways of construct- 
ing this sort of object more explicitly (see Refs. [23,37,38] for example), but here it is useful 
to exploit the construction used in the classical case. The amplitude to pass through A is 
therefore given by 

^a(x/, t|xo, 0) = y Vx{t) exp {iS[x{t)]) 9 dt /A(x(t)) - (3.23) 

Here, the ^-function ensures that only paths x{t) that spend a time in excess of e in A 
contribute to the sum. 

The class operator Ca is a propagator at fixed energy, E, say, so this is given by 

/oo 
rfr e-*^^^a(x/,T|xo,0) (3.24) 
-oo 

When ga is replaced with an unrestricted propagator, we require that (3.24) is annihilated 
hj H — E, and this is why we choose an infinite range for r, rather than a half-infinite 
one [39,40,41]. (As we shall see below, when ga is a restricted propagator, we encounter 
some difficulties here, although the correct range for r is still the infinite one). The total 
(regularized) decoherence functional is therefore given by 

//■oo roo 
d"a;/d"a;oc^"4 / dr dr' e''^^ e'^ ^ 
.1 — so V — oo 

X ^a(x/,T|xo,0) ^*,(x;,t'|xo,0) *£;o(xo)*|, (x[,) (3.25) 
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This may also be written, 

/oo rOO 
dr / dr' e-'^^ e'^^ 
-oo J —oo 

X J Vx{t) J Vx'{t)exp(^iS^[x{t)]-iS^'[x'{t)]^ 

X ^Eoi^oWE'^o) (3.26) 

where note that the two actions in the path integral are over different ranges of time. It is 
straightforward to show that 

D{a, a') = 5{E - Eq)5{E' - E'q)5{Eq - E'^) (3.27) 

In the induced inner product scheme we therefore replace the righthand side by 1, verifying 
that the construction is correctly normalized. 

3(D). Modified Class Operators 

The basic scheme described above runs into a interesting difficulty in that the class 
operators defined by (3.24) do not satisfy the constraint equation. We have, for example, 
for the class operator for paths that enter the region A, 

Ca(x/,xo) = P dre-'^^ J Vx{t) exp {iS[x{t)]) 9 dt fA{^{t)) - (3.28) 

It may be shown that this satisfies the constraint everywhere except on the boundaries of 
the region A. That is, there is a discontinuity as one of the end-points crosses the boundary 
of A. This is an issue because a basic rule of the game of constructing the decoherence 
functional for these systems is that we only work with objects which satisfy the constraint 
equation (or operators which commute with it). Mathematically, this is to ensure that the 
underlying symmetry, reparametrization invariance, is fully respected (the induced inner 
product, note, is defined only between objects which satisfy the constraint). Physically, 
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it is related to the fact that the universe is a closed system, and measurements of it (the 
class operators are generalizations of the notions of measurement) must not displace its 
wave function. 

Because of this difficulty, it is necessary to define a modified class operator C'^ which is 
as close as possible to the path integral one, but satisfies the constraint equation everywhere 
[22] . A way to do this is to first compute the class operator when one of the end-points 
is inside A. This defines a solution to the constraints for one end-point inside and the other 
outside. We then define C'^ to be the object which satisfies the constraint everywhere and 
matches this expression for one end-point inside and the other outside. 

An example of this difficulty was encountered and solved in Ref . [23] , which concerned 
the decoherent histories analysis of the relativistic particle. Suppose one considers the 
class operator for propagating between spacetime points Xq and x^p with the restriction 
that the paths never cross the spacelike surface x^ = constant. This object, denoted 
Cr{xj,XQ), is readily constructed using the free propagator together with the method of 
images. However, it does not satisfy the constraint everywhere: the constraint operating 
on Cr gives a (5-function on the spacelike surface. One can also see physically why there 
might be problems here. Suppose for simplicity that Xj and xq are timelike separated 
and consider the trajectory of a classical particle passing through these two points. Such 
a trajectory, if extended beyond these points, must cross every surface of constant x^ . 
Similarly, in the quantum theory, there are no non-trivial solutions to the Klein- Gordon 
equation that are zero on one side of a spacelike surface. Hence the only sensible possible 
answer for the class operator for this situation is Cr = 0. This is indeed the solution used 
in Ref. [23] and led to physically expected results. 

This issue concerning the replacement of the original class operators with modified ones 
is related to the fact that in the original class operator (3.28), the functions restricting 
the paths to enter the region A involve a time integral over a finite interval [0, r], whereas 
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the classical result (2.13) involves a similar restricting function but with a time integral 
over an infinite range. If a finite range is used in the classical result, it is no longer 
reparametrization invariant, as we saw in Section 2. This is why we discussed the classical 
result at such length. 

The construction of these modified class operators is crucial to the construction of the 
decoherence functional for reparametrization-invariant theories. It is probably fair to say 
that the method suggested above for constructing them has not at this stage been fully 
explored. We will show below that a physically plausible modified class operator is readily 
constructed in the semiclassical approximation, but a more thorough investigation of this 
issue will be deferred to a later publication. It should also be noted that reparametrization- 
invariance can be rather subtle. For example, the path-integral constructed class opera- 
tors Eq.(3.28) certainly appear to be reparametrization-invariant at the level of symmetry 
transformations at the Lagrangian level, yet do not quite satisfy the constraint. Here, we 
have used the expression "reparametrization-invariant" to mean satisfying the constraint 
everywhere (or having zero Poisson bracket with the Hamiltonian everywhere, in the clas- 
sical case). This issue is related to the connections between Lagrangian and Hamiltonian 
symmetries, and between the path integral and Dirac quantizations. See Refs. [21,41] for 
further discussion. 



4. THE SEMICLASSICAL APPROXIMATION TO THE 
DECOHERENCE FUNCTIONAL AND PROBABILITIES 

We may now compute the decoherence functional. It is 

Dia,a') = J d^'xfcPxocPyo C^(x;,xo)C2(x;,yo) *(xo)**(yo) (4.1) 

(For convenience, here and in what follows we will in fact drop the notation involving 

23 



different values of E to regularize the expressions, unless necessary). To see how the 
semiclassical approximation works out, we will in this Section assume decoherence (for 
example, by restricting to special initial states) and concentrate on the construction of the 
probabilities that the system passes through the region A. We will return to the question 
of decoherence for general initial states in Section 6. 

We begin by computing the semiclassical approximation to the modified class operator 
C^(xy, xq). Recall that it is given by a suitable modification of the path integral expression 
(3.28) with one end-point inside A and the other outside. In the absence of any restrictions 
on the paths, the path integral will be dominated by the classical paths connecting the 
initial and final points. The classical paths will be the solutions to the equations of motion 



Mi + VV(x) = 



(4.2) 



which satisfy the boundary conditions 



x(0) = XQ, x(t) = X/ 



(4.3) 



In addition, these paths must satisfy the constraint equation 




(4.4) 



This equation determines the time r in terms of xq, Xj- and E, hence the final form of the 
extremizing paths have no reference to time. It is also useful to introduce ^(x^,xo), the 
classical action from xg to xj. It obeys the time- independent Hamilton- Jacobi equation 



with respect to each end point. 



1 



{VAf + V = E 



(4.5) 



2M 



The initial and final momenta are given by derivatives of the classical action. 



Pf = '^fM'^f^^o)^ PO = -Vo^(x;,xo) 



(4.6) 
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The semiclassical approximation to the unrestricted path integral is given by a sum of 
terms each of the form 

G(xy,xo)=P(x^,xo)e^^("/'"o) (4.7) 

The quantity P is a prefactor, whose specific form is lengthy to calculate, but will not in 
fact be required. (For the case of the time-dependent propagator, the prefactor would be 
given in terms of the determinant of the matrix of second derivatives of ^(xj-,xo). Here, 
because of the constraint equation (4.5), this matrix is in fact singular and the expression 
for the prefactor is more complicated [42,43]). 

The semiclassical form (4.7) satisfies the constraint equation (i.e., is annihilated by 
H — E) in the WKB approximation, as it should. If the classical trajectory from xq to x^- 
is unique there will be just one term in the semiclassical approximation. If there is more 
than one, there will be a sum of similar terms, one for each trajectory. For the moment 
we will assume that there is just one. We will also assume that the extremizing classical 
solution is real, and thus the action ^(xj-,xo) is real. 

With the restriction that the paths must pass through A, we expect that the class 
operator will be given again by (4.7) when the classical path passes through A, and will 
be zero when the classical path does not pass through A. It is then not difficult to see 
that the modified class operator for this case may therefore be written, 

C^(x;,xo) = {^j'^Jt fA{4(t)) - Pi^f^^) e^^(^/'^°) (4.8) 

The ^-function here is the same as in the (rewritten) classical case Eq.(2.13) in terms of 
"initial" and "final" points, where x^(t) denotes that classical path from xq to x^-. (This 
is exactly as in the classical case depicted in Fig.l). Note also that 

• V ^(^P dt /A(x^(t)) - = (4.9) 

as may be shown by shifting the t integration. It follows that the modified class operator 
is a WKB solution to the constraint equation, as required. We have therefore succeeded 
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in computing, in the semiclassical approximation, a modified class operator satisfying the 
constraint equation everywhere, corresponding to the restriction to paths passing through 
the region A. This is a very simple result but turns out to be crucial to the rest of the 
derivation. 

It is important that t is integrated over an infinite range in the quantity inside the 
^-function, otherwise the modified class operator would not in fact satisfy the constraint. 
Recall that the originally defined class operator Eq.(3.24) contained a similar ^-function, 
with a finite range of time integration, which one might have been tempted to use in the 
semiclassical approximation, but this class operator does not in fact satisfy the constraint. 

Hence we see that the difference between the modified and original class operators in the 
semiclassical approximation is the difference between using the entire classical trajectory 
or using finite segments of it in the 6'-functions. We also see that these modified class 
operators are the correct ones to use in order to be consistent with the discussion of the 
classical case and Eq.(2.13). There, we saw that it is appropriate to sum over classical 
paths intersecting A even if A does not lie on the segment of classical trajectory between 
xq and Xf. This feature therefore appears to be necessary for the particular type of 
reparametrization invariance used here. Only the entire trajectory is reparametrization- 
invariant notion. A finite section of trajectory is not. (See Ref.[41] for a further discussion 
of reparametrizations in this sort of context). 

The off-diagonal terms of the decoherence functional are now given in the semiclassical 
approximation by 



It is now essentially a sum over pairs of classical paths, and the ^-functions restrict the 




X P(xy,xo)P*(xj,yo) exp(iA(xj,xo) -z^(xj,yo)) p(xo,yo) (4.10) 
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paths to either pass on not pass through the region A. We see also from the invariance 
property of the ^-functions, (4.9), that they are invariant under shifting the regions Aq, 
along their classical trajectories. This is the expression of the idea that the decoherence 
functional, in the semiclassical approximation, knows only about entire trajectories. 

It is now convenient to introduce the variables 

1 

Xo=2(^0 + yo), v = xo-yo (4.11) 



and thus 



xo = Xo + ^v, yo = Xo - (4.12) 



and we also rewrite the density operator in terms of the Wigner function Eq.(3.11). We will 
discuss the detailed mechanism of decoherence in the next Section. For the moment, we 
will simply assume decoherence, which essentially means assuming that v is concentrated 
around zero, and work out the form of the probabilities. Although we note that this 
assumption can be justified for initial states p(xo, yo) which are approximately diagonal 
in position. We now set v = in the prefactors P and in the ^-functions, and we obtain 
the probability 

PA = J d'^Xfci'^Xorvdy (^y"°° dt /A(X^(t)) - \P{yif,Xof 

X exp (^a(xy,Xo + ^v)-a(x^,Xo-^v)+ip-v^ W^(p,Xo) (4.13) 

Expanding the action terms to linear order in v (decoherence again allows us to drop the 
higher order terms), the v integral may be performed and we obtain 

PA = j d^Xfd'^Xod^p dt /A(X^(t)) - 

X |P(x^,Xo)|^5(") (p + VoA(xy,Xo)) W^(p,Xo) (4.14) 

where Vq operates on the initial point Xq. Finally, the integration over xy may be per- 
formed. The delta-function constraint then means that the quantity X^(t) (the classical 
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path from Xq to Xj-) is replaced by X^'(t) (the classical path with initial data Xo,po)- 
Although we have not worked out the explicit form of the prefactor P, we deduce that it 
must in fact drop out when the Xj^ integration is carried out, because the probability must 
equal 1 when the ^-function is removed. (From this we deduce that {Pi'x.j, Xo)P must be 
the Jacobian factor in the change of integration variables from xj to Vo^(xy,Xo)). We 
therefore obtain the final result, 

PA = j d^'Xody e (^J^^ dt /A(X^^(t)) - Wip, Xo) (4.15) 

As expected, this is the classical result Eq.(2.10) with the classical phase space distribution 
function replaced by the Wigner function. 



5. SYSTEMS OF HARMONIC OSCILLATORS 

There is one simple system for which the discussion of decoherence and probabilities 
is particularly simple, and this is the case of a collection of harmonic oscillators. It also 
enjoys the property that its spectrum is discrete, hence the induced inner product is not 
required for normalization. 

In non-relativistic quantum mechanics, in the search for emergent quasi-classical his- 
tories, it is of interest to consider histories characterized by strings of phase space quasi- 
projectors Pp. These are positive hermitian operators concentrated on a region F of phase 
space, but are not quite projectors since they only have Pp ?a Pp. Omnes has proved an im- 
portant theorem about these projectors, which is essentially that they are approximately 
preserved in form under unitary evolution and moreover approximately follow classical 
evolution [12]. That is, 

e-'^^Pre'^^ ^ Pr, (5.1) 
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where Tt is the classical evolution of the phase space cell F. The approximation holds 
when the phase space cells are significantly larger than a quantum sized cell, and for times 
not so long that wave packet spreading becomes significant. In the special case of the 
harmonic oscillator, the approximation holds for all time. This result allows one to show 
that firstly, histories of phase space projectors are approximately decoherent for a wide 
variety of initial states, and secondly, that their probabilities are peaked about classical 
evolution. Differently put, on sufficiently coarse grained scales, quantum systems have an 
approximate determinism that ensures decoherence and approximate correspondence with 
classical physics. 

It seems reasonable to suppose that the timeless models considered here might have 
analogous properties. We will demonstrate this for a system of harmonic oscillators in an 
energy eigenstate. The Hamiltonian for a set of N identical harmonic oscillators is 



factors). The standard coherent states (see Ref.[44], for example) are denoted |p,x) and 
they have the important property that they are preserved in form under unitary evolution. 



where Pt,^t cire the classical solutions matching p,x at t = 0, hence they are strongly 
peaked about the classical path. In Ref.[27], a set of states were introduced which are 
timeless analogues of the usual coherent states. They are 




(5.2) 



(the factor of is included to subtract the vacuum state energy and avoid certain phase 




(5.3) 



\ct>p^) = 6{Ho-E)\p,^) 




(5.4) 
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These states were referred to in Ref.[27] as "timeless coherent states". (See also [14,45]). 
They are exact eigenstates of the Hamiltonian, 



Furthermore, since the coherent states |pi,Xi) are concentrated at a phase space point for 
each t, integrating t over a whole period produces a state which is concentrated along the 
entire classical trajectory. They are therefore the natural analogues of the usual coherent 
states. Their properties are similar in many ways to the usual coherent states and are 
described in more detail in Ref.[27]. 

Each state is labeled by a fiducial phase space point p, x which determines the classical 
trajectory the state is peaked about. Under evolution of the fiducial point p, x to another 
point, ps, Xg, say, along the same classical trajectory, the state changes by a phase. 



Two timeless coherent states of different energy are exactly orthogonal. If they have the 
same energy then they are approximately orthogonal if they correspond to sufficiently 
distinct classical solutions. They also obey a completeness relation. 



[Note that the notation S{Hq — E) is a rather loose one. This object is really the projection 
operator onto the subspace of energy E, for which it is exactly true that S{Hq — E)'^ = 
5{Hq — E)]. Since 5{Hq — E)\'ip) = lip) on any solution to the eigenvalue equation {Hq — 
E) = 0, this is essentially a completeness relation on the set of solutions to the eigenvalue 
equation. We may therefore write any solution |^) as a superposition of timeless coherent 
states. 



(5.5) 



|0px)-|</>p.x.) = e^^^|</>px) 



(5.6) 




(5.7) 




(5.8) 
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Given these preliminaries, we may now discuss the decoherence functional. We will 
consider coarse grainings in which the paths in configuration space either pass or do not 
pass through a series of regions denoted A = Ai, A2 • • •, we will take Ai, A2 • • • to lie along 
a classical path. Hence we need at least two such regions to fix a configuration space path. 

The decoherence functional, in terms of the modified class operators, is, in the semi- 
classical approximation, 

i;(A, A) = J d^xf^xoryo C^(x;,xo)Cg(x;,yo) *(xo)**(yo) (5.9) 

The modified class operator C^(xj-, xq) is given by Eq.(4.8), so is equal to the unrestricted 
semiclassical propagator G when xj or xq lie on the classical path specified by A and is 
zero otherwise. Also, = 5{Hq — E) — C'^. 

We first consider the case in which the initial state is a timeless coherent state, \4>px)- 
It is then straightforward to see that 

c;^i<^px) ~ i<^px) (5.10) 

when the trajectory labeled by the fiducial points p, x passes through the regions A, and 

C'^IM^O (5.11) 

otherwise. Also, since 

5(iyo-£^)|0px) = |0px) (5.12) 

it follows that 

C^|</'px)~0 (5.13) 

when the trajectory labeled by p, x passes through A. From these results it is easy to see 
that the decoherence functional is approximately diagonal. Furthermore, the probability 
for entering the regions A is then approximately 1 or 0, depending on whether the classical 
trajectory of the timeless coherent state passes through A. 
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Now consider the case of a more general initial state. We expand it in timeless coherent 
states, as in Eq.(5.8). Using the above results, we therefore find 

r jN jN 

C'^\i^) ^ l</'px)(0px|V) (5.14) 

Here, D denotes the set of phase space points p, x whose classical trajectories pass through 
the regions A in configuration space. Similarly, 

CaIV') « ^ ^§1^ l0px)(0px|V) (5.15) 

where D denotes the set of phase space points whose classical trajectories never pass 
through A. Clearly we again have approximate decoherence, because of the approximate 
determinism. We may therefore assign a probability for passing through A, 

It is the integral over the phase space region A of the phase space distribution function 
[(^pxlV')! • Because is an eigenstate of the Hamiltonian, it is easy to see using the 
definition (5.4) of the timeless coherent states that 

('/'pxlV') = (px|V) (5.17) 

and so the probability now is 

PA - / ^^^^|(px|V^)P (5.18) 
JD (27r)^^ 

It is then a standard result that the integrand is in fact a smeared Wigner function 

|(px|V)|2 = j d^p'd^x' e-i(P-P')^-|(x-x')^ I^(p',x') (5.19) 

This object is positive even though the original Wigner function W of is not [44]. 
Hence we obtain a result which is essentially identical to the classically anticipated result 
(2.10), with a smeared Wigner function as the phase space distribution function. 
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It is also of interest to note that the result for the probability may be written in the 
form 

PA = {^\Pd\^) (5.20) 
where Pd is the approximate projection operator 

r jN jN 

Moreover, since the timeless coherent states |0px) are exact eigenstates of i^o, we have 
that 

[Pd.Hq]=Q (5.22) 

Hence, we see that the result may be written in the standard quantum-mechanical form 
for a probability, in terms of an operator which commutes with the constraint. 

The result here, of approximate decoherence and simple expressions for the probabil- 
ities, like the corresponding non-relativistic result is due to the approximate determinism 
contained in the quantum theory. It works only when we ask for the probabilities for 
approximately classical histories. To obtain probabilities for more complicated histories, 
and for systems which are not harmonic oscillators (where there is wave packet spreading), 
we need an environment to produce decoherence. 



6. DECOHERENCE THROUGH AN ENVIRONMENT 

As stated above, the decoherence functional is typically not diagonal for most initial 
states, and a physical mechanism is required to produced decoherence. In this section we 
therefore consider the addition of an environment to produce decoherence of histories. The 
results of this Section therefore simply justify the assumed decoherence of Section 4, and 
little affect the final result of the probabilities, but it is important to see in detail how this 
works. 
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6(A). Semiclassical Approximation to the 
Decoherence Functional with Environment 



For what we will do here, the specific form of the environment turns out not to be 
very important. But for definiteness, we take the environment to be a large collection of 
harmonic oscillators with coordinates denoted q^^, where A runs over a large number of 
values, with a linear coupling to the system. For notational simplicity we will assume that 
for each system coordinate x in the n-dimensional configuration space there is a set of n 
oscillators with coordinate q for the environment. The case of of more oscillators is easily 
obtained from this. The total action of the system is 

5 = 5o[x]+5,[x,q] (6.1) 

and the corresponding Hamiltonian is 

iy = i?0(x) + iy,(x,q) (6.2) 
We shall assume that the state of the whole system has the form 

*(x,q)=V(x)x(x,q) (6.3) 

This may be inserted into the Wheeler-DeWitt equation {H — E)^ = to obtain a pertur- 
bative solution about the solution with no environment. We will concentrate on the case 
in which the wave function is of oscillatory form, so the background solution is of WKB 
form 

V'(x) = C(x)e*'^W (6.4) 
where S obeys the Hamilton-Jacobi equation 

^ {VSf + V(x) = E (6.5) 
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and C obeys the equation 

V^C + 2VS' • VC = (6.6) 
The environment wave functions obey the Schrodinger equation 

iVS ■Vx = H,x (6.7) 

We will consider the case of a superposition of WKB states in Section 7. 
The decoherence functional now is 



"OO , , 



/oc roo 
dr / dr 
-00 ^ — 00 

X / Px(t) / Vy{t)exp(iS^[^{t)]-iS^'[y{t)] 
Ja Jo! ^ 

X F[x(t),y(t),T,T'] V£;o('^o)V'|^(yo) (6-8) 

(suspending for the moment the necessity to use modified class operators). Here 5'q denotes 
the action over the fixed time range [0, r] (and note that the time ranges are different on 
either side of the decoherence functional). The influence functional F[x(t), y(t), r, r'] is 
given by, 

F[x(t),y(t),T,T'] = y" Pq(t)Pr(t)exp(z5|[x,q] -z5^'[y,r]) 

X x(xo,qo)x*(yo,ro) (6.9) 

It is different in form in two ways to the usual influence functional. Firstly, the time ranges 
on either side are not the same, and secondly, the initial state x depends on the system 
variables x, with a different dependence on either side of the influence functional. The 
functional integral is over all pairs of paths q(t), r(t) which meet at the flnal point 

q(T) = r(T') (6.10) 

and this point is summed over. The paths also match the initial values qq, tq, which are 
then folded into the initial state. 
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It is useful to go now to the semiclassical approximation for the system variables. We 
also now recall that we must use modified class operators for the system variables (this 
does not affect the environment dynamics, at this level of approximation). The off-diagonal 
terms of the decoherence functional are now given by 

D{A, A) = J d"xjd"xod"yo (^J^ dt /A(xi(t)) - 

X ^'(x/,xo)P*(x^,yo) exp (z^(xy,xo) -z^(x^,yo)) 
X F(x_^,xo,yo) V'(xo)V'*(yo) (6.11) 
Here, F{^f, xq, yo) is the influence functional with the semiclassical approximation for the 
system variables inserted. That is, for x(t) we insert the classical trajectory from xg to xj 
in time r, and the value of r is then determined (in terms of xq and xy) by the constraint 
equation, and similarly for y{t). The decoherence functional is again essentially a sum 
over pairs of classical paths for the system variables, the path from xq to xj, and the path 
from yo to xy. 

6(B). Calculation of the Influence Functional 

We may now calculate the influence functional with the semiclassical approximation 
for system variables inserted. The influence functional may be written, 

^(x/,xo,yo) = y" d^'g/ 0(q/,xj,xo)0*(q/,x_^,yo) (6.12) 

where 

(/•(q/, xj, xq) = j d"?o f (q/, x/|qo, xq) x(qo, xq) (6.i3) 

Here, we have introduced the propagator for the environment variables along the system 
classical trajectory from xg to xy, 

^(qj,xj|qo,xo) = J Vq{t) exp {iS}[^{t),q{t)]) (6.14) 
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The environment state x ni^iy be normalized according to 

J rq |x(q,x)|2 = l (6.15) 

for all X. Since g propagates x unitarily along a fixed system trajectory, it follows that 

J d'^qf |0(q/,xj,xo|^ = l (6.16) 

This means that the influence functional satisfies 

|F(xj,xo,yo)|^<l (6.17) 

with equality when xq = yo, indicating that the influence functional is peaked about 
xq = yo, which is the decoherence effect we need. 

The influence functional is difficult to evaluate in general. It can be evaluated exactly 
when both g and x Gaussian. It is also effectively Gaussian in form if we have a large 
number of oscillators in the environment, for then the many-oscillator inffuence functional 
is essentially the original one, (6.9), raised to a high power. This strongly enhances the 
peaking about xg = yo- A simple and reasonably general form for the inffuence functional 
may therefore be obtained by expanding about xo = yo, and assuming either a Gaussian 
form, or a large number of oscillators (or both). 

We again use the coordinates Xo, v deffned in Eq.(4.11). We have 

1 1 
0(q/,x/,xo) = (^(q/,xy,Xo) + -v" aa(^(q_^,xj.,Xo) + -v%^dadf,(j){cif,yif,'KQ) + ■■■ 

(6.18) 

where da — d/dX^. Inserting in the inffuence functional, and also introducing the notation 
00 = 0(q/,x/,Xo), we get 

+ \v''v^ j rqf {(P^QdadbcPo + cPodadf^^Q - 2dact>Qd^,ct>l) 

+ ••• (6.19) 
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By differentiating the normalization of 0, (6.16), we see that 

j (Pqf <^Ml = - j (Tqf daMl (6-20) 

and 

j rqf {<^Qdadk<^l + <^ldadM = ' j d^'qf (^a^o^ft^o + 56<^o5a0^) (6.21) 
Using these relations, we may now write the influence functional as 

^ exp (^iv^Ta - \v^v^ {E^t, - r^r^,)^ + • • • (6.22) 
The coeflBcients Ta and E^^^ are given by 

ra(xy,Xo) ='-J (f'qf {^oda^l - ^Ida^o) 

S^;,(xy, Xo) = j dTqf da(t>Ml (6-23) 

The approximation of writing F as a Gaussian becomes exact when g and x are Gaus- 
sian, and is also true approximately when there are a large number of oscillators in the 
environment. We have therefore obtained the influence functional, as required. 

6(C). Reparametrization Invariance 
in the Influence Functional 

The form of the influence functional indicates, as expected, that there is a suppression 
of interference for paths with xq / yo- This is the usual decoherence effect. However, the 
situation in the reparametrization-invariant theory considered here is not so simple. We 
expect that the decoherence functional depends, in some sense, only on reparametrization- 
invariant quantities, and in the semiclassical approximation used here, this means it de- 
pends on entire classical paths (rather than individual points). Differently put, we do not 
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expect (or need) the destruction of interference for points xq, yo lying on the same classical 
path (connecting each to xy), since these points are effectively equivalent. What we expect 
is that the influence functional will not be exponentially small when xy, xq, yo he along a 
single classical path. We must therefore see how reparametrization invariance is expressed 
in the influence functional. 

We first note that 

00 = y" d'^qo £f(q/,xy|qo,Xo) x(qo,Xo) (6.24) 

and let us see how this quantity varies with Xq. If this were the usual non-relativistic 
quantum mechanics, with propagation from initial time to to final time tj, then (pQ would 
in fact be independent of tQ. We expect a similar property here, that is, that (f)Q is constant 
(as a function of Xq) along a certain vector field. The initial state x(qO)Xo) obeys the 
Schrodinger equation 

iWoS{Xo) ■ Vox = i?£(qo,x/,Xo)x (6.25) 

The propagator g on the other hand, obeys Schrodinger equations with respect to both 
the final and initial points, 

zV/^(x/,Xo) • Vfg = H,{qf,yif,Xo)g (6.26) 
iVoA(xj,Xo) • Vq^ = H,{qQ,^f,Xo)g (6.27) 

(Note that the expected minus sign in the Schrodinger equation with respect to the initial 
point is already contained through the fact that Vo^ is minus the initial momentum, as in 
Eq.(4.6)). Now the point is here that g and x obey different Schrodinger equations, so at 
this stage, 0o does not obviously have any constant directions in Xo - neither Vo'S' • Vq^o 
nor Vo^ • Vo0o zero. 

However, as we saw in the Section 4 (without environment), the path integral enforces 
the condition p = — Vq^. We anticipate that this condition is approximately enforced 
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with the environment in place. Furthermore, the initial Wigner function for a WKB state 
is of the approximate form, 

l^(p,Xo) = |C(Xo)|2 5(p- Vo^(Xo)) (6.28) 

It follows that the sum over paths is dominated by configurations for which 

VoA(xj,Xo)^-Vo^ (6.29) 

This means that the trajectories of S are the same as the classical trajectories from Xq to 
xj. From this it is then easy to show that 

Vo^(Xo) • V(/)o ~ (6.30) 

(essentially for the same reason that the analogous non-relativistic version is independent 
of the initial time to)- I^i the influence functional, two neighbouring points xq, yo on the 
same classical trajectory have v = xq — yo proportional to Vq-S'. It follows that 

v'^Ta = 0, v'^J^ab = (6.31) 

which means that the influence functional does not suppress interference between points 
on the same trajectories, only between points on different trajectories. That is, when the 
condition p = — Vq^ is true, we get the expected result that the influence functional is a 
function only of entire trajectories, and not of the individual points along those trajectories. 

6(D). Decoherence and the Evaluation of the v integral 

The off-diagonal terms of the decoherence functional may now be written 
D{A,A) = J (Pxf(PXQ(Pvdy 9 (^J^ dt /A(x^(t)) - 

\J—oo / 
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xP(x;,Xo + ^v)P*(x;,Xo-^v) 

X exp(i(VoA(x;,Xo)+p)-v + 0(v3)) l^(p,Xo) 

X exp (^iv'^Ta - lv%^a^,^ (6.32) 

where 

(^ab = ^ab " ^a^b (6.33) 

and we have again introduced the variables Xq, v in the exponential part. The decoherence 
functional is a sum over pairs of classical paths, one set of paths intersecting xo,Xj- and 
passing through A at any stage along the path, the other set of paths intersecting yo,xy 
and never passing through A at any stage along the path (see Fig.l). It is easily seen that 
for this particular coarse graining, in which we are interested in paths that either pass or 
do not pass through the region A, we do not in fact encounter the situation discussed in 
Section 6(C), in which xg, yo, xj lie along the same classical trajectory. That is, the coarse 
graining is such that v is never proportional to Vq^, and the potentially singular situation 
(6.31) does not arise. The influence functional therefore does its job of suppressing the 
contribution from non-zero values of v to a decoherence width determined by the inverse of 
the non-zero eigenvalues of cTq^. If the size of the coarse graining region A is greater than 
this width then the off-diagonal terms of the decoherence functional are approximately 
zero. 

We therefore have approximate decoherence and we may examine the probability for 
passing through A, which is 

p{A) = J (rxf(rXQ(rv(rp dt /A(X^(t)) - 

xP(x^,Xo + iv)P*(x^,Xo-^v) 
X exp(z(Vo^(xy,Xo)+p)-v + 0(v3)) I^(p,Xo) 
X exp (^iv^Ta - ^v^v^aab^ (6.34) 
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where v has been set to zero in the ^-function. We are now summing over pairs of classical 
paths which both pass through the region A, so now we do have the possibility of xq, yo, xy 
lying along the same path, and hence the matrix u^j, is potentially singular, by (6.31). This 
means that some care is necessary in the v integral. 

If we formally carry out the integral over v, we get, 

PA=J rxfrxody e dt /A(x^(t)) - 

X |P(x;,Xo)|V(p,Xo) 

X exp (^-^(VoA(x;,Xo) + p + r)V-HvoA(x^,Xo) + p + r)^) (6.35) 

Changing integration variables from xj to pg = — Vo^(xj^,Xo) This is conveniently writ- 
ten 

PA = y d^Qd^XQ e (^£2 dt /A(X^^(t)) - W^(po, Xq) (6.36) 
where we have defined the smeared Wigner function 

w^(po,Xo) = j dy exp(^-^(po-p-r)V-i(po-p-r)^^ w-(p,Xo) (6.37) 

Again the prefactors P drop out in the change of integration variables, as in Section 4. This 
smearing of the Wigner function represents environmentally induced fluctuations about the 
classical evolution, and the small additional term F in the exponent represents the back 
reaction of the environment on the classical equations of motion. 

As stated, the results (6.31) suggest that the matrix a^i, is singular, but it is easy to 
see the significance of this. When the matrix is non-singular, the v integral produces a 
Gaussian peak about po = P + T, which, as we have seen represents fluctuations about 
classical evolution. If the matrix is singular in a certain direction, it is easy to see from 
(6.34) that the v integral in this direction will produce a (5-function, instead of a Gaussian. 
It will still be peaked about the same conflguration, but there are no fluctuations in that 
direction. 
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The result for the probabihties therefore approximately coincides with the classical 
result Eq.(2.10). We have concentrated on the case in which A is a single region of 
configuration space, but the result straightforwardly generalizes to the case in which A 
consists of a series of regions Ai, A2, • • •. The above result then shows that the probability 
is peaked when the series of regions lies along a classical path (plus the environmental 
effects of a small back reaction and small fluctuations). 



7. SUPERPOSITION STATES 

The calculations of Sections 4 and 6 concerned only single WKB states of the form 
(6.4). It is therefore important to reconsider the decoherence calculation of Section 6 for 
the more general case of a superposition of WKB states, 

= $1 + $2 = Cie'^'xi + C2e'^'x2 (7.1) 

This turns out in fact to be quite straightforward, mainly because similar calculations 
(involving the reduced density matrix, not the decoherence functional) have already been 
done. 

Inserting Eq.(7.1) in the decoherence functional, we obtain, a result of the form 

D = Dii + D12 + D21 + D22 (7.2) 

where, in an obvious notation, Dn is the decoherence functional with initial density matrix 
|\E'l)(\E'l|, D12 is the decoherence functional with the operator |\E'i)(\E'2| in the initial state 
slot, and so on. Clearly the analysis of Dn and D22 is identical to the case considered 
already - we get decoherence, and probabilities given in terms of the Wigner functions of 
^'l and '^2- Hence these two terms correspond to a statistical mixture of the two initial 
states. 
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The interesting terms are D12 and D21 (= 0^2)^ which correspond to interferences 
between different WKB branches. Prom Section 6, we see that 

/OO roo 
-00 J— 00 

X I V^{t) J Vy{t)e^p{iS^[^{t)]-iS^'[y{t)]) 

X Fi2[x(0,y(t),T,T'] Ci(xo)e^^i("o) C2*(yo)e-^^2(yo) (7.3) 

where 

i^l2[x(t),y(t),T,T'] = J Pq(t)Pr(t)exp (i5^[x, q] - z5|'[y, r]) 

X xi(xo,qo)x2(yo,ro) (7.4) 

(As in Section 6(A) we should then replace the class operators with their modified version 
and go to the semiclassical limit). Now it is easy to see that the influence functional is 
the overlap of the two initial states, but with each unitarily evolved along two different 
trajectories. That is, in the semiclassical approximation for the system, 

F12 = (X2(yo)|C/'''(x/,yo)t/(x/,xo)|xi(xo)) (7.5) 

where t^(xj, xq) denotes the unitary evolution of the environment states along the system 
classical trajectory from xq to xy. 

Clearly |-Fi2|^ < 1, and because F12 is an overlap between a pair of states it will 
typically be such that F12 is strictly less than 1. In this case, when raised to a high power, 
as happens when we take a large number of oscillators in the environment, we will get a 
very strong suppression of terms with xq 7^ yg- In particular, even when xg = yo, we get 

i^l2 = (X2(xo)|xi(xo)) (7.6) 

which will be less than 1, quite simply because xi and X2 are different states. We therefore 
find that D12 and D21 are much smaller than the diagonal terms Dn, D22- This destruc- 
tion of interference between WKB states therefore comes about for essentially the same 
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reason that the corresponding off-diagonal terms of the density matrix are very small, as 
discussed previously [46,47,48,49]. 

It should be noted, that in Eq.(7.5), it could in fact happen that F12 = 1, as a result 
of a careful choice of xi 7^ X2 together with a suitable choice of xg 7^ yo, in particular, if 

t/(x/,xo)|xi(xo)) = C^(x/,yo)|x2(yo)) (7.7) 

The point here, however, is that this becomes very unlikely with a large environment. 
With a large collection of oscillators in the environment, the environment states are a 
tensor product over A of states |x^(xo)), for example, and then 

F12 = n (x^(yo)|t^^(x/,yo)C/(x/,xo)|x^(xo)) (7.8) 
A 

As the size of the environment goes to infinity, the possibility of Eq.(7.8) being exactly 1 
becomes negligible. 

It is also of interest to look at the special case in which the wave function Eq.(7.1) 
is real (as is the case in the no-boundary wave function of Hartle and Hawking), so that 
*2 = a^nd X2 = Xi- When xq = yo, we then have that 
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\Fi2r 



(7.9) 



dqo xi(xo,qo) 

Since xi is generally complex (it obeys the complex Schrodinger equation (6.7)), the right- 
hand side of Eq.(7.9) will again be less than 1, so the argument still goes through [49]. The 
argument fails if xi is real. But then it would have to be an eigenstate of the environment 
Hamiltonian for all values of xq, and this would not lead to decoherence, so we may 
disregard this case. 



8. SUMMARY AND DISCUSSION 



We have studied the quantization and interpretation of simple timeless models de- 
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scribed by an equation of the type (1.1). In particular, we studied the question, what is 
the probabihty that the system passes through a region A of configuration space without 
reference to time? 

We obtained the classical answer to this problem, in three different forms, in terms of a 
classical phase space distribution function t(;(p,x), satisfying {H,w} = 0, the analogue of 
(1.1). This function needs some care in normalization since it is constant along the (possibly 
infinite) classical orbits. A very useful step in the classical case was the introduction of the 
phase space quantity Eq.(2.7), which is (5-function peaked on the classical path and also 
has vanishing Poisson bracket with the Hamiltonian, so is an observable. This quantity 
assists in understanding some aspects of the quantum theory. 

We constructed the decoherence functional, following the general scheme of Refs. [21,22], 
using the induced inner product. Although the general scheme has been presented pre- 
viously, a key part of our contribution to this area is the explicit identification (in the 
semiclassical approximation) of the class operator Eq.(4.8) satisfying the constraint every- 
where describing histories restricted to pass through a region A of configuration space. 
Having made this identification, a major part of our work was to show that the decoherent 
histories approach then reduces, approximately, to the corresponding classical result, but 
with the classical phase space distribution function w{p,:x.) replaced by the Wigner func- 
tion VF(p, x) of the quantum theory. We also explored the decoherence and probabilities 
for a system of harmonic oscillators using the timeless coherent states, in terms of which 
the analysis is particularly transparent and fully agrees with intuitive expectations. 

In brief, therefore, we have shown that heuristic classically inspired notions of inter- 
pretation for simple timeless models may in fact be derived from the decoherent histories 
analysis of such models. This result is by no means unexpected, but the key aspects of the 
derivation are the elucidation of the role of the constraint and the related reparametriza- 
tion invariance in the construction of both the classical and quantum results. Furthermore, 
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the complete absence of a time parameter is not an obstruction to quantization. 

There are a number of issues which the present work has generated and will be discussed 
in a later publication, but we mention them briefly below. 

First of all, the main difficulty in computing the decoherence functional for our chosen 
coarse graining is the calculation of the modified class operators. Even before modification 
class operators of the type Eq.(3.21) are difficult to calculate (typically they can only be 
obtained exactly in the very simple situations where the method of images may be used). 
The suggested scheme for constructing modified operators obeying the constraint has not 
yet been explored fully. Here, we have constructed physically plausible modified class 
operators in the semiclassical approximation, obtaining full agreement with the classical 
results. Some exact modified class operators for simple coarse grainings of the relativistic 
particle have been constructed in Ref. [23] , but it is not yet clear how general those results 
are. Hence a more detailed investigation of these modified class operators is called for. 
(We note in passing that, from the simple examples in which the modified class operators 
have been calculated, their calculation does in fact appear to be considerably easier than 
the original ones, Eq.(3.21)). 

Second, we have assumed (except in Section 5) that both the initial states and the prop- 
agators are in the oscillatory regime. This means in the propagator that we assume the 
dominant contribution comes from real configurations (rather than Euclidean or complex 
ones). Many interesting models in quantum cosmology have a Euclidean region, corre- 
sponding, for example, to "tunneling from nothing". It is not immediately clear how the 
semiclassical calculation of Section 4 is modified to include this case, the main difficulty 
being understanding what the class operators are. This case is therefore probably related 
to the question of a more general formula for the modified class operators. 

Third, it is generally understood that decoherence of histories is related to the existence 
of "records" [10,33]. This means that it is possible to find a projection operator Ra which 

47 



is perfectly correlated with the class operators Ca in terms of which the probabilities may- 
be written, 

p{a) = Tr (CapCa) = TV (Rap) (8.1) 

In the case of a non-relativistic model where decoherence is produced by an environment, 
it is possible to explicitly identify the environmental variables which store the information 
about the system [33]. It would be very desirable to do this in the timeless case consid- 
ered here. It seems likely that the variables are very similar to the case of Ref.[33], but 
the interesting question is the role of reparametrization invariance in this situation, and 
whether the records are closely related to observables in the operator approach. 

A fourth issue concerns the connection between the decoherent histories analysis con- 
sidered here and the master equation for the reduced density operator (in the case where 
decoherence is produced by an environment). It would be of interest to see if the discussion 
of decoherence and probabilities can be re-expressed in the simpler language of the density 
operator, as it sometimes can in non-relativistic decoherence models. This is currently 
under investigation [50]. 

These and related issues will be taken up in future publications. 
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Figure 1 

The rewritten classical probability Eq.(2.13) in terms of a sum over initial and final 
points xq and x^-. The probability for not entering A is a sum over paths as in case (a). 
The probability for entering A includes a sum over classical paths in which A lies between 
the initial and final points, as in case (b). But, to agree with the phase space form of the 
result Eq.(2.10), it must also include a sum over initial and final points for which A does 
not lie between them, as in case (c). This figure also applies to the semiclassical propagator 
Eq.(4.8) in the quantum case. 
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